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As a continuation of a systematic study of reactions relevant to the astrophysical p process, the 
cross sections of the 74,76 Se(p,7) 75 ' 77 Br and 82 Se(p,n) 82 Br reactions have been measured at energies 
from 1.3 to 3.6 MeV using an activation technique. The results are compared to the predictions of 
Hauser-Feshbach statistical model calculations using the NON-SMOKER and MOST codes. The 
sensitivity of the calculations to variations in the optical proton potential and the nuclear level 
density was studied. Good agreement between theoretical and experimental reaction rates was 
found for the reactions 74 Se(p,7) 75 Br and 82 Se(p,n) 82 Br. 



PACS numbers: 25.40.Lw, 26.30.+k, 27.50.+e, 97.10.Tk 



I. INTRODUCTION 



Details of the nature of the astrophysical p process 
0, @ producing the rare, proton-rich, stable isotopes of 
elements beyond Ni still remain elusive. It has been 
shown that such proton-rich isotopes can be produced 
by photodisintegrations in massive stars 0, Q , involving 
(7,n), (7,p), and (7, a) reactions at stellar temperatures 
exceeding 10 9 K. Depending on intricate details of the 
stellar structure and evolution, p nuclides are mostly pro- 
duced in the final explosion of a massive star (M > 8M Q ) 
as a core-collapse supernova or in pre-explosive oxygen 
burning episodes 0,3- Despite the fact that p nuclei 
can be produced consistently with solar ratios over a wide 
range of nuclei in such scenarios, there remain deficien- 
cies in a few regions, for mass numbers A < 124 and 
168 < A < 200. The problem is especially pronounced 
in the Mo-Ru region where the p isotopes are strongly 
underproduced. This fact motivates the search for alter- 
native or additional ways to produce these nuclides. Pro- 
ton capture and photodisintegration processes in explod- 
ing carbon-oxygen white dwarfs have been suggested as a 
source 0, El j or thermonuclear explosions in the proton- 
rich layer accreted on the surface of a neutron star in 
a binary system with mass flow from the main-sequence 
companion star 0, • 

Common to all approaches is that the modelling re- 
quires a large network of hundreds of nuclear reactions 
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involving stable nuclei as well as unstable, proton-rich 
nuclides. It is well known that (7,0:) reactions are im- 
portant at higher masses, whereas (7,p) or proton cap- 
ture is important for the production of the less massive 
p nuclei. The relevant astrophysical reaction rates cal- 
culated from the cross sections are inputs to this net- 
work, therefore their knowledge is essential for p process 
calculations. While there are compilations of neutron 
capture data along the line of stability, there are still 
very few charged-particle cross sections determined ex- 
perimentally llllSlIllIlEQHEini despite big 
experimental efforts in recent years. Thus, the p process 
rates involving charged projectiles are still based mainly 
on (largely untested) theoretical cross sections obtained 
from Hauser-Feshbach statistical model calculations. 

The aim of our systematic study is to contribute to 
the existing database of measured cross sections relevant 
to the astrophysical p process, and to check the reliabil- 
ity of the statistical model calculations over an extensive 
set of nuclides. This way the uncertainties in the p pro- 
cess abundance calculations arising from nuclear physics 
input can be constrained. 

In this paper we present measurements of the 
74 Se(p, 7 ) 75 Br, 76 Se(p, 7 ) 77 Br, and 82 Se(p,n) 82 Br reac- 
tions in the astrophysically relevant energy range using 
an activation technique. The two proton capture reac- 
tions and their inverses are directly playing a role in the 
synthesis of 74 Se, whereas the (p,n) reaction can be used 
as a further test of statistical model calculations. The 
choice of these reactions is further elaborated in Sec. |n] 
and the experimental method is described in Sec. II I II 
The resulting cross sections and astrophysical S factors 
(given in Sec. IIV|I are compared with results obtained 
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with the two Hauser-Feshbach statistical model codes 
NON-SMOKER [3 13 HI and MOST [2l| in Sec. |V| 
Both codes are based on the statistical theory of Hauser 
and Feshbach [^ but use different models for the nuclear 
ingredients of the calculations. By using both codes we 
are able to investigate the effects of a broader range of 
nuclear level densities and optical model potentials. The 
astrophysical reaction rates derived from our new exper- 
imental data are given in the concluding Sec. IVII 



II. INVESTIGATED REACTIONS 

The element Se has six stable isotopes with mass num- 
bers A = 74, 76, 77, 78, 80 and 82 having isotopic 
abundances of 0.89%, 9.36%, 7.63%, 23.78%, 49.61% 
and 8.73%, respectively. Proton capture reactions on 
these isotopes lead to Br isotopes among which 79 Br 
and 81 Br are stable. Therefore the cross sections of 
the 78 ' 80 Se(p,7) 79 ' 81 Br reactions cannot be measured us- 
ing an activation technique. The half life of 78 Br, 
i.e. the reaction product of 77 Se(p,7) 78 Br, is too short 
(T1/2 = 6.49 min.) for our experimental method (see 
experimental details) . Thus the aim of the measurement 
was to determine the proton capture cross section of three 
reactions: 74,76,82 Se(p,7) 75 ' 77,83 Br. However, in the case 
of 82 Se the (p,n) channel opens already at E\ a b — 891 keV 
and, therefore, competes strongly with the (p,7) reac- 
tion in the whole energy range investigated in the present 
work. Moreover, the strongest 7 transition following the 
(3 decay of 83 Br has a very low relative intensity (only 
1.2 % of the decays lead to the emission of this 529.6 keV 
7-ray). The latter two facts made it impossible to observe 
the decay of 83 Br and to measure the 82 Se(p,7) 83 Br cross 
section. However, the 82 Se(p,n) 82 Br cross section could 
be determined. 

In summary, the cross sections of three reactions have 
been measured: 74 ' 76 Se(p,7) 75J7 Br, and 82 Se(p,n) 82 Br. 
The astrophysically relevant energy range (Gamow- 
window) for these reactions in the temperature range 
from T = 1.8xl0 9 K to 3.3xl0 9 K spans from 1.25 to 
3.87 MeV. This energy region was covered by the exper- 
iment. The relevant part of the chart of nuclides can be 
seen in Fig. ^ where the proton induced reactions and 
the decay of the reaction products can also be seen. The 
decay parameters used for the analysis are summarized 
in Table □ 
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FIG. 1: The relevant part of the chart of nuclides with the 
decay of the reaction products. The stable isotopes are rep- 
resented by bold squares. 



TABLE I: Decay parameters of the Br product nuclei taken 
from literature. 



Product 
nucleus 


Half life 
[hour] 


Gamma 
energy 
[keV] 


Relative 
intensity 
per decay 
\%] 


Reference 


/b Br 


1.612 ± 0.002 


286.6 
141.3 


88 ± 5 
6.6 ± 0.5 


w 


77 Br 


57.036 ± 0.006 


239.0 
520.6 
297.2 


23.1 ± 0.5 
22.4 ± 0.6 
4.16 ± 0.21 


[24J 


82 Br 


35.3 ± 0.02 


776.5 
554.3 
619.1 


83.5 ± 0.8 
70.8 ± 0.7 
43.4 ± 0.4 


[25] 



activation procedure. Aluminum is ideal as backing ma- 
terial because no long-lived activity is produced during 
its bombardment with protons in the investigated energy 
range. Moreover, Al can easily be distinguished from the 
Se in the RBS spectrum (RBS was used to monitor the 
target stability, see below). 

The target thickness was measured with the proton 
induced X-ray emission (PIXE) technique at the PIXE 
set-up of the ATOMKI According to PIXE re- 

sults, the target thickness was ranging from 200 to 
700/ig/cm 2 , corresponding to a proton energy loss of 
lOkeV (at 3.6MeV) to 60keV (at 1.3MeV), respectively. 
The thicker targets (500-700 /ig/cm 2 ) were used at the 
lower part of the bombarding energy range, where the low 
cross section results in very low induced 7-activity. The 
proton energy loss has been calculated with the SRIM 
code Wi- 



lli. EXPERIMENTAL PROCEDURE 

A. Target properties 

The targets were made by evaporating metallic Se with 
natural isotopic abundance onto a thick Al backing. Nat- 
ural targets have the advantage that the three investi- 
gated reactions can be studied simultaneously in a single 



B. Activation 

The activations were carried out at the 5 MV Van 
de Graaff accelerator of the ATOMKI by irradiating the 
Se targets with a proton beam. The energy range from 
E p =1.3 to 3.6 MeV was covered with 100 - 300 keV steps. 
The schematic view of the target chamber can be seen in 
Fig. |21 After the last beam defining aperture the whole 
chamber served as a Faraday-cup to collect the accumu- 
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FIG. 2: Schematic view of the target chamber. 



lated charge. A secondary electron suppression voltage of 
—300 V was applied at the entrance of the chamber. Each 
irradiation lasted about 10 hours with a beam current of 
typically 5 to 10 /iA. Thus, the collected charge varied 
between 180 and 360 mC. The current was kept as stable 
as possible but to follow the changes the current integra- 
tor counts were recorded in multichannel scaling mode, 
stepping the channel in every minute. This recorded cur- 
rent integrator spectrum was then used for the analysis 
solving the differential equation of the population and 
decay of the reaction products numerically (see eqs. 5-10 
inrefB3). 

A surface barrier detector was built into the chamber 
at 9=150° relative to the beam direction to detect the 
backscattered protons and to monitor the target stability 
this way. The RBS spectra were taken continuously and 
stored regularly during the irradiation. In those cases 
when target deterioration was found, the irradiation was 
repeated with another target. Final cross section results 
were derived only from those measurements where no 
noticable target deterioration was found, i. e. the growth 
of the Se peak area in the RBS spectrum was uniform 
with collected charge within the statistical uncertainty 
(below 1%). 

The beam was wobbled across the last diaphragm to 
have a uniformly irradiated spot of diameter of 8 mm on 
the target. The target backing was directly water cooled 
with an isolated water circulating system. 

Between the irradiation and 7-counting, a waiting time 
of 30 minutes was inserted in order to let the disturbing 
short lived activities decay. 



C. Detection of induced 7-radiation 

The 7 radiation following the /3-decay of the produced 
Br isotopes was measured with a HPGe detector of 40% 
relative efficiency. The target was mounted in a holder at 
a distance of 10 cm from the end of the detector cap. The 
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FIG. 3: Activation 7 spectrum after irradiating a target with 
2.5 MeV protons. The 7 lines used for the analysis are in- 
dicated by arrows. The not labeled peaks are either from 
laboratory background or beam-induced background on im- 
purities in the target and backing or 7 lines from Se + p 
reaction products which were not used for the analysis. 



whole system was shielded by 10 cm thick lead against 
laboratory background. 

The 7 spectra were taken for at least 12 hours and 
stored regularly in order to follow the decay of the dif- 
ferent reaction products. 

The absolute efficiency of the detector was measured 
with calibrated 133 Ba, 60 Co and 152 Eu sources in the 
same geometry used for the measurement. Effect of the 
finite sample size (beam spot of 8 mm in diameter) was 
measured by moving the point-like calibration sources 
over this surface and measuring the difference in effi- 
ciency. This was then included in the 7% error of detector 
efficiency. 

Fig. 13 shows an off-line 7-spectrum taken after irra- 
diation with 2.5 MeV protons in the first lh counting 
interval. The 7 lines used for the analysis are indicated 
by arrows. 

Taking into account the detector efficiency and the rel- 
ative intensity of the emitted gamma-rays, coincidence 
summing effects were for all three reactions below 1% 
and were neglected. 



IV. EXPERIMENTAL RESULTS 

Tables ITT1 - IIVI summarize the experimental results for 
the three studied reactions. The quoted errors of the 
center-of-mass energies given in the first column corre- 
spond to the energy loss in the targets calculated with 
the SRIM code. The error of the cross section (S factor) 
values is the quadratic sum of the following partial errors: 
efficiency of the HPGe detector (~7%), number of target 
atoms (~6%), current measurement (3%), uncertainty of 
the level parameters found in literature (<6%), counting 
statistics (0.1 to 15%). 

At E p — 2. 176 MeV proton bombarding energy the 
77 Se(p,n) 77 Br reaction channel opens, having the same 
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TABLE II: Experimental cross section and S factor of the 
74 Se(p,7) 75 Br reaction 



TABLE IV: Experimental cross section and 5* factor of the 
82 Se(p,n) 82 Br reaction 



E c 



Cross section 



[keV] 




[/ibarn] 


1 A £ K 
1400 


± 


27 


0.0 


± 


U. ^ 


1009 


i 

± 


12 


1 A Q 

1U.S 


i 

± 


1.9 


loOo 


i 

± 


22 


on q 


i 

± 


o.z 


1766 


± 


12 


41.7 


± 


6.7 


1858 


± 


zu 


60.7 


± 


9.0 


1954 


± 


22 


94.0 


± 


14 


2057 


± 


18 


119 


± 


18 


2159 


± 


14 


191 


± 


28 


2310 


± 


11 


366 


± 


55 


2456 


± 


13 


467 


± 


69 


2752 


± 


13 


1039 


± 


154 


3049 


± 


12 


1561 


± 


231 


3348 


± 


8 


1888 


± 


282 


3547 


± 


6 


2966 


± 


483 



S factor 
[10 6 MeV-b] 



5.6 
6.6 
6.4 
6.1 
5.0 
4.4 
3.2 
3.1 
3.0 
2.1 
1.6 
0.98 
0.54 
0.53 



1.1 

1.2 

1.0 

1.0 

0.8 

0.7 

0.5 

0.5 

0.5 

0.3 

0.2 

0.15 

0.08 

0.09 



E c , 


m. 




Cross section 


5* factor 


[keV] 




[/ibarn] 


no 6 


MeV-b] 


IzOo 


zh 


31 


U. 10 


i 

± 


U.14 


8.2 


i 

± 


i c: 
1 .0 


140 ( 


i 

zh 


27 


0.0 


i 

zh 


U.o 


5.6 


± 


U.o 


10 1 1 


zh 


12 


1 1 Q 

11.8 


i 

± 


1 9 

1.2 


7.2 


zh 


U. / 


lOOU 


zh 


22 


1 7 
1 / .0 


i 

± 


1. ( 


5.5 


i 

zh 


U.o 


1768 


± 


12 


44.1 


± 


4.3 


6.5 


± 


0.6 


1860 


± 


20 


79.8 


± 


7.6 


6.6 


zb 


0.6 


1956 


± 


22 


137 


± 


13 


6.5 


± 


0.6 


2060 


± 




219 


± 


21 


5.9 


± 


0.6 


2162 


± 


14 


367 


± 


35 


6.0 


zb 


0.6 


2313 


± 


11 


742 


± 


71 


6.1 


zb 


0.6 


2459 


± 


13 


1142 


± 


109 


5.1 


zb 


0.5 


2755 


± 


13 


3330 


± 


317 


5.1 


± 


0.5 


3052 


± 


12 


6778 


± 


647 


4.2 


± 


0.4 


3352 


± 


8 


12448 


± 


1208 


3.6 


± 


0.4 


3552 


± 


6 


19329 


± 


2257 


3.5 


± 


0.4 



TABLE III: Experimental cross section and S factor of the 
76 Se(p,7) 77 Br reaction. Note that above 2.148 MeV cm. 
energy the 77 Se(p,n) 77 Br reaction has a contribution to the 
quoted experimental values. See text for details. 



Ec.m. 

[keV] 




Cross section 
[/ibarn] 


S factor 
[10 6 MeV-b] 


1456 


± 


27 


9.5 


zb 


0.8 


14.8 


zb 1.3 


1569 


± 


12 


10.4 


± 


0.9 


6.3 


=b 0.6 


1658 


± 


22 


21.8 


zb 


2.2 


6.8 


=b 0.7 


1767 


± 


12 


59.6 


zb 


6.2 


8.8 


=b 0.9 


1859 


± 


20 


53.8 


± 


5.3 


4.4 


zb 0.4 


1954 


± 


22 


93.1 


zb 


9.2 


4.4 


zb 0.4 


2058 


± 


18 


147 


± 


14 


4.0 


zb 0.4 


2160 


± 


14 


207 


zb 


20 


3.4 


=b 0.3 


2311 


± 


11 


796 


zb 


78 


6.5 


=b 0.7 


2457 


± 


13 


1182 


zb 


116 


5.3 


=b 0.5 


2753 


± 


13 


3200 


zb 


312 


4.9 


=b 0.5 


3050 


zb 


12 


6467 


± 


631 


4.0 


zb 0.4 


3349 


zb 


8 


11529 


zb 


1143 


3.3 


=b 0.3 


3548 


± 


6 


17252 


zb 


2044 


3.1 


zb 0.4 



final nucleus as 76 Se(p,7) 77 Br. Thus, above this energy 
these two reactions cannot be distinguished with the ac- 
tivation technique using targets with natural isotopic 
abundance. The measured cross section and S-factor 
values above the threshold are actually the weighted 



sum of the two cross sections: a cxp = <J\ + 0.82 -<72, where 
(7 1 and (72 are the cross sections of 76 Se(p,7) 77 Br and 
77 Se(p,n) 77 Br, respectively and 0.82 stands for the iso- 
topic ratio of 77 Se and 76 Se: 7.63%/9.36%. 



V. COMPARISON WITH PREDICTIONS 



each figure. Two different combinations of nuclear level 
densities (NLD) and optical potentials (OMP) are used 
by each code, the OMP of |2l with the NLD of |1 
in MOST, and the OMP of HI with the NLD of H3 
in NON-SMOKER. From the figures it is apparent that 
for the Se isotopes studied herein, the results of the for- 
mer combination are systematically higher (by a factor 
of 2) than the latter at all studied energies. Concerning 
the energy dependence, a decrease at low energy is seen 
for 82 Se(p,n) 82 Br which seems not to be present in the 
data. The predictions of the latter input combination 
are in very good agreement with the data for energies 
larger than 1.7 MeV in 74 Se(p,7) 75 Br, but cannot repro- 
duce the observed energy dependence at smaller energies. 
They are higher than experiment by factors 1.25-2 for 
76 Se(p,7) 77 Br above 1.8 MeV and cannot reproduce the 
experimental variation in the S factor at lower energies. 
They are able to reproduce the data very well for the 
82 Se(p,n) 82 Br reaction above 1.7 MeV. Again, the fea- 
tures in the experimental S factors below that energy 
cannot be reproduced. However, the deviations do not 
exceed a factor of 1.5. As will become evident in Sec. lVIl 
these deviations at the lower end of the measured energy 
range do not contribute significantly to the astrophysical 
reaction rate for the temperatures given there. 

In summary, contrary to some previous measurements 
of (p,7) reactions (see e.g. [l5 | L the standard predictions 
using the nuclear inputs of |3Ct l3l| are in good agree- 
ment with the present experimental data and most of 
other proton capture data 1 . Even more so, if one con- 
siders that the calculations use global parameters which 



In Figs. 01-0 we show the comparison of the experi- 
mental data with theoretical results. At first, we discuss 
the standard predictions of the Hauser-Feshbach codes 
MOST and NON-SMOKER, shown in the part (a) of 



Note: Ref. ll'J points out that the theoretical values in I 111 were 
mistakenly given too high. Using the co rrec t value, there is good 
agreement also with the experiment of 1 1 (I . 
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are not fine-tuned to the specific nuclei involved. For 
such global calculations, a 30% deviation averaged over 
all nuclei is not unusual, similarly to what was found for 
neutron capture |3lj . Locally, larger errors are possible, 
of course. Within this range of uncertainty the predic- 
tions agree very well, especially regarding the astrophys- 
ical rates as given in Sec. lVIl However, as can already be 
seen from the comparison of the standard predictions, the 
results are sensitive to the nuclear inputs, even though 
the reaction mechanism, and thus the reaction model to 
be applied, is unambiguous. Therefore, to better under- 
stand the contributions of different model inputs it is 
necessary to disentangle the effects of mainly the NLD 
and the OMP. This will also help to constrain future, 
improved parametrizations of these nuclear properties. 
For this purpose, additional calculations were performed 
using both codes and varying the NLD and OMP mod- 
els. The combination of NLDs and OMPs used in each 
calculation are summarized in Table IVl In the following 
sections, we discuss each reaction separately. 



A. 74 Se(p, 7 ) 75 Br 

The comparison of the results for different NLDs and 
OMPs is shown in the parts (b) and (c) of Fig. 0] The 
labels are explained in Table Ivl It should be noted that 
INP-1 and INP-5 are the default predictions also shown 
on part (a) in Fig.H Moreover, INP-4 and INP-10 both 
use the same level density and optical potential. The re- 
maining difference between them has to be attributed to 
further differences in other nuclear inputs (such as the 
photon width) because the basic approach is the same in 
both calculations. One has to keep in mind this addi- 
tional small difference when comparing results obtained 
from the two codes (results INP-1 to INP-4 and INP-5 
to INP-10, respectively). 

As can be seen from the figure, the dependence on the 
OMP is much stronger than the sensitivity to the NLD, 
especially at the lower end of the studied energy range. 
All of the calculations with the same OMP but differ- 
ent NLDs show a similar energy dependence whereas a 
change in the potential not only leads to different abso- 
lute values but also to different excitation functions. At 
the low energies studied here, with no other open chan- 
nels than (p/y), the cross section depends only on the 
incident proton transmission coefficients, the NLD in the 
compound nucleus determining the photon widths, and 
the NLD in the target entering the compound-elastic pro- 
ton channel. The effect of the NLD in the latter channel 
is small, as transitions with higher proton energies to low- 
lying states will be dominating and because these states 
are experimentally and explicitly included in the calcu- 
lation. From Fig. 0] part (c) it is evident that none of 
the OMPs are able to reproduce the data at low ener- 
gies, except for the ones of 32] (INP-3) and [p3 (INP- 
6). The equivalent square- well potential of |34j | and the 
phenomenological OMP of [12 are able to give the most 
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FIG. 4: Experimental S factor of the reaction 74 Se(p,7) 75 Br 
and statistical model calculations. The standard predictions 
of MOST (INP-1) and NON-SMOKER (INP-5) are compared 
to the data (a). Below, the effect of varyingthe nuclear level 
density with fixed proton optical potential |30l| is shown (b), 
as well as the effect of varying the optical potential (c). In 
the part (c) of the figure, the curves labeled INP-1 to INP-3 
include the level density of [2s|, the ones labeled INP-5 to 
INP-8 the one of 31]. See Table IVl for an explanation of the 
labels. 
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TABLE V: Overview of the inputs in the Hauser-Feshbach calculations. 



Label 


Level density 




Optical potential 


INP-1 " " 


Demetriou & Goriely (2001) 


[29] 


Bauge et al. (2001) [28] 


INP-2 " 


Demetriou & Goriely (2001) 


[29] 


Jeukenne et al. (1977) [301 


INP-3 a 


Demetriou & Goriely (2001) 


[29] 


Koning (2002) [32] 


INP-4 " 


Thielemann et al. (1986) [33] 




Jeukenne et al. (1977) [301 


INP-5 c d 


Rauscher et al. (1997) [31] 




Jeukenne et al. (1977) [30] 


INP-6 c 


Rauscher et al. (1997) (311 




equivalent square well [34] 


INP-7 c 


Rauscher et al. (1997) (311 




Becchetti & Greenlees (1969) [351 


INP-8 c 


Rauscher et al. (1997) [3JJ 




Perey (1963) [£6J 


INP-9 c 


Holmes et al. (1976) [341 




Jeukenne et al. (1977) [301 


INP-10 c 


Thielemann et al. (1986) \^} 




Jeukenne et al. (1977) [£Q] 



""Calculated with MOST 

'Standard prediction as available at 

http: //www- astro .ulb. ac .be/Html/hfr .html (version 
09/12/2002) 

Calculated with NON-SMOKER 

^Standard prediction as published in llit |20H : also available at 
: //nucastro . org/reaclib . html 



satisfactory overall description of the data. The latter 
phenomenological potential has been obtained for ener- 
gies from 1 keV up to 200 MeV, using a unique, flat 
functional form for the energy dependence. This may be 
the reason why it gives a relatively flatter variation with 
energy compared to the other OMPs. What is somewhat 
surprising, is the good agreement obtained with the sim- 
ple equivalent square-well potential. 

The OMPs with the largest deviation are the ones of 
(INP-7) and [H (INP-1). The former OMP has been 
derived for energies above 20 MeV, so it is not surprising 
that it fails. The latter is based on microscopic nuclear 
matter calculations to start with, but has been subse- 
quently re-adjusted to reproduce an extensive database 
of reaction observables. However, the reliability of the re- 
sulting renormalization coefficients has only been tested 
for energies in the 10-30 MeV range. The overpredic- 
tion of the data could thus be attributed to inappropriate 
values of the renormalization coefficients at low energies. 
Another important effect might be caused by deforma- 
tion which we further discuss in Sec. IVII 



B. 76 Se(p, 7 ) 77 Br 

In this case, the data cannot be straightforwardly com- 
pared with the predictions because of the complication 
with the 77 Se(p,n) channel. As discussed in Sec. II VI the 
experiment was not able to distinguish between 76 Se(p,7) 
and 77 Se(p,n) because they have the same final nucleus 
77 Br. Therefore, only the points measured below the 
(p,n) threshold at £^ = 2.176 MeV, are included in the 
comparison. 

The effect of varying the NLDs and the OMPs is shown 
in Fig. [S] The labels are explained in Table IVl What was 
said for 74 Se(p,7) applies equally here. Again, the depen- 
dence on the OMP is the strongest one and a large spread 



in absolute values and energy dependence is found. The 
OMPs of [35| and [2^| significantly overestimate the data, 
for the possible reasons mentioned in the previous sec- 
tion. The OMPs of 34] and give a very flat energy 
dependence which describes the data well over most of 
the energy region. On the other hand, the slightly steeper 
excitation functions obtained with the OMPs of [3j| and 
[30| are able to give an overall reasonable account of the 
data. 

The low-energy structure found in the experiment (first 
and fourth data point at the lowest energies) cannot 
be reproduced by any of the OMPs and it is doubtful 
whether such a behavior of the S factor could be found 
in a statistical model calculation at all. 



C. 82 Se(p,n) 82 Br 

The effect of different NLDs and OMPs is shown in 
Fig. HJ The labels are explained in Table What was 
said for 74 Se(p,7) in Sec. IV Al applies similarly here. 

For this case, the dependence on the optical potential 
is much stronger than the sensitivity to the NLDs over 
the whole energy region. The neutron emission channel 
opens at the energy of 0.8 MeV and rapidly becomes the 
most dominant channel at all measured energies. In such 
a case, the HF cross section depends mainly on the in- 
cident proton transmission coefficients which is exactly 
what is observed in the figures. Similar to the results ob- 
tained for the other Se isotopes, the OMP of [H (INP-1) 
and [3j| (INP-7) overpredict the data by at least a fac- 
tor of 2. Here, also the potential of [3|| (INP-8) yields 
a significantly higher S factor although the shape of the 
energy dependence agrees well. Much better is the agree- 
ment of the OMP of [3^] (INP-3). However, its energy 
dependence is slightly too flat. Again, very good agree- 
ment with the data over the whole energy range is found 
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with the equivalent square-well potential of |3J| (INP- 
6). The microscopic potential of [H (INP-2, INP-5) de- 
scribes well the data above 1.7 MeV but shows a different 
energy dependence below that energy. 
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FIG. 5: Experimental S factor of the reaction 76 Se(p,7) 77 Br 
and the statistical model calculations. The standard predic- 
tions of MOST (INP-1) and NON-SMOKER (INP-5) are com- 
pared to the data (a). Below, the effect of varying the nuclear 
level density with fixed proton optical potential |3(Jl is shown 
(b), as well as the effect of varying the optical potential (c). 
In the part (c) of the figure, the curves labeled INP-1 to INP- 
3 include the level density of [2sJ, the ones labeled INP-5 to 
INP-8 the one of [3l|. See Table IVI for an explanation of the 
labels. 



VI. SUMMARY AND CONCLUSIONS 

We have measured three proton-induced reac- 
tions in the astrophysically relevant energy range: 
74 ' 76 Se(p,7) 75 ' 77 Br, and 82 Se(p,n) 82 Br. Using an activa- 
tion technique we were able to obtain reaction cross sec- 
tions and S factors in the energy range relevant to the as- 
trophysical p process. The new data were compared with 
predictions of the Hauser-Feshbach statistical theory us- 
ing global models of OMPs and NLDs. An extensive 
investigation of the sensitivity of the theoretical calcula- 
tions to the different inputs was presented. As expected, 
in this astrophysically relevant low-energy region the re- 
sults are more sensitive to the OMP than to the NLD. 
The best overall agreement is obtained with the global 
microscopic optical potential of [3(j and an equivalent 
square- well potential j34[. The latter, somewhat surpris- 
ing, observation seems to be in contradiction with a 
vious comparison of rates using these two potentials 
where the largest deviations between them were found in 
the mass range A > 70. However, this is strongly energy 
dependent since in Fig. 4 of |3jj the largest deviations 
arise at even lower proton energies and the differences 
vanish quickly at higher energies. Likewise, we expect to 
find larger differences between the predictions with the 
two potentials when going to even lower energy. 

Further good agreement with absolute values of the 
present data is found with the potential of [12 • How- 
ever, the obtained energy dependence of the S factor re- 
mains slightly too flat. This is caused by the functional 
dependence on the energy assumed in the definition of 
this potential. 

The Se isotopes considered here are deformed and 
therefore deformation effects might be important. For 
the calculations using the OMPs of [30l l35| the well- 
known approach of employing a spherical equivalent to a 
deformed potential (with a larger diffuseness parameter) 
was used. However, from our analyses above it can be 
seen that even a purely spherical potential, such as the 
equivalent square well, can yield satisfactory results. On 
the other hand, deformation is important for interpreting 
the results obtained with the OMP of [2£| ■ A significant 
fraction of the shown deviations is due to the spherical 
treatment. A further investigation of dependencies on 
the OMP parameters and how to include deformation ef- 
fects in OMPs, and more particularly coupled channels 
calculations, is beyond the scope of this paper. 

Reaction rates are the relevant quantities in astrophys- 
ical applications. In Tables IVIHVIIII we give such as- 
trophysical reaction rates computed from our data in- 
cluding errors. The temperature range for each reac- 
tion was chosen numerically so that all significant contri- 
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TABLE VI: Experimental reaction rates computed from the 
S factor of the 74 Se(p,7) 75 Br reaction and comparison to pre- 
dicted values. 



T 9 
[10 9 K] 


Experimental rate 
[cm 3 s" 1 mol -1 ] 


Prediction [201 
[cm 3 s -1 mol - 


2.00 


2.720X10 1 


± 


5.301 


3.75X10 1 


2.50 


2.177xl0 2 


± 


3.978X10 1 


2.55xl0 2 


3.00 


9.454 xlO 2 


± 


1.663xl0 2 


1.02xl0 3 


3.50 


2.813xl0 3 


± 


4.843 xlO 2 


2.92 xlO 3 


4.00 


6.490xl0 3 


± 


1.104xl0 3 


6.75xl0 3 


4.50 


1.251 xlO 4 


± 


2.115xl0 3 


1.33 xlO 4 


5.00 


2.117xl0 4 


± 


3.566xl0 3 


2.36 xlO 4 


6.00 


4.616xl0 4 


± 


7. 764 xlO 3 


5.79 xlO 4 




(b) 



2000 2400 

Ec.m.tkeV] 




2 - 

(C) 

o I 1 

1200 1600 



2000 2400 2800 

Ec.m. [keV] 



FIG. 6: Experimental S factor of the reaction Se(p,n) Br 
and the statistical model calculations. The standard predic- 
tions of MOST (INP-1) and NON-SMOKER (INP-5) are com- 
pared to the data (a). Below, the effect of varying the nuclear 
level density with fixed proton optical potential |3Cll is shown 
(b), as well as the effect of varying the optical potential (c). 
In the part (c) of the figure, the curves labeled INP-1 to INP- 
3 include the level density of [2S|, the ones labeled INP-5 to 
INP-8 the one of [3j]]. See Table IVI for an explanation of the 
labels. 



TABLE VII: Experimental reaction rates computed from the 
S factor of the reactions 76 Se(p,7) 77 Br and 77 Se(p,n) 77 Br, and 
comparison to predicted values. 



Tg 


Experimental rate 


Prediction [20J 






°"(p,7) + 0.82(T( P:n ) 


[10 9 K] 


[cm 3 s -1 mol -1 ] 


[cm 3 s" 1 mol" 1 ] 


2.50 


4.540x10^ ± 5.671 xlO 1 


6.05 xlO 2 


3.00 


2.189xl0 3 ±2.590xl0 2 


3.17x10 s 


3.50 


7.237x10 s ±8.343xl0 2 


1.14xl0 4 


4.00 


1.821 xlO 4 ± 2.078x10 s 


3.18xl0 4 


4.50 


3.758xl0 4 ±4.275xl0 3 


7.42 xlO 4 



butions to the integration over the proton energy came 
from within the energy range of our data. Also shown 
is a comparison to the rates from the standard NON- 
SMOKER prediction [2(j; reaction rates for the cases 
INP-1 to INP-10 can be derived considering the fact 
that rates from slowly varying S factors scale approxi- 
mately with the S factor. Excellent agreement is found 
for 74 Se(p,7) 75 Br and 82 Se(p,n) 82 Br. This illustrates 
how deviations with respect to the data are averaged 
out by the integration involved in the calculation of the 
reaction rate, especially at the edge of the respective 
Gamow window. The prediction overestimates the rate 
of 76 Se(p,7) 77 Br + 77 Se(p,n) 77 Br ( Table IVTTt by factors 
of 1.18-1.77, with better agreement at lower tempera- 
ture. Thus, the present work confirms the trend seen 
in previous investigations of proton-induced reactions for 
intermediate mass targets, namely that there is overall 
acceptable or good agreement between data and global 
predictions. Apart from a few cases where some devi- 
ations seem to persist independent of nuclear input (Sr 
isotopes studied in |l5j), the discrepancies between the- 
oretical calculations are not as large as those observed 
for different optical a potentials (see, e.g., [TTL EM l39| ). 
This seems to hold for the mass range 70 < A < 100 
HHIHHE0E1EII3- Nevertheless, more 
measurements are required to completely cover the rele- 
vant mass region and provide constraints on the nuclear 
input used in Hauser-Feshbach calculations. 
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TABLE VIII: Experimental reaction rates computed from the 
S factor of the 82 Se(p,n) 82 Br reaction and comparison to pre- 
dicted values. 



Tg 

[10 9 K] 


Experimental rate 
[cm 3 s -1 mol -1 ] 


Prediction [20] 
[cm 3 s -1 moP] 


1.50 


1.610 


± 1.541xl0 _i 


2.07 xl0 u 


2.00 


4.561X10 1 


± 4.639 


4.91 xlO 1 


2.50 


4.377xl0 2 


± 4.515X10 1 


4.45 xlO 2 


3.00 


2.231xl0 3 


± 2.333xl0 2 


2.32 xlO 3 


3.50 


7.519xl0 3 


± 7.982 xlO 2 


8.52xl0 3 


4.00 


1.907xl0 4 


± 2.054xl0 3 


2.45 xlO 4 
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